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ABSTRACT 

We propose a novel approach to the investigation of how black hole spin influences the behavior of astrophysical accretion. 
We examine the terminal behavior of general relativistic, multi-transonic, hydrodynamic advective accretion in the Kerr metric, 
by analyzing various dynamical as well as thermodynamic properties of accreting matter extremely close to the event horizon, as 
a function of Kerr parameter a. The examination of the transonic properties of such flows can be useful in self-consistently dis- 
criminating between prograde and retrograde relativistic accretion, and in studying the spectral signatures of black hole rotation. 
Subject headings: accretion, accretion disks — black hole physics — hydrodynamics — relativity 



1. INTRODUCTION 

General relativistic (GR) accretion flows dive supersonic ally 
onto astrophysical black holes (BH) in order to satisfy the in- 
ner boundary condition imposed by the event horizon. Hence 
such accretion flows must be globally transonic (e.g., Fukue 
1987; Chakrabarti 1989). This extremely hot, high density, 
ultra-fast flow can provide key features of the spectra of some 
galactic and extragalactic BH candidates. It is believed that 
most of the supermassive BHs powering active galactic nuclei 
(AGN) and stellar mass BHs in galactic microquasars possess 
non-zero values of the BH spin parameter, a, which might 
play a crucial role in determining their spectra and jet pro- 
duction (e.g., Rees 1984; Frank, King & Raine 1992; Kato, 
Fukue & Mineshige 1998). Hence it is very important to study 
the dynamical as well as the thermodynamic properties of GR 
transonic BH accretion very close to the event horizon and to 
investigate how BH spin modulates such behavior 

Study of the multi-transonic behavior of BH accretion was 
initiated by Abramowicz & Zurek (1981). Here we formulate 
and solve the conservation equations governing GR, multi- 
transonic, advective accretion in the Kerr metric to calculate 
the flow behavior extremely close (~ 0.01 gravitational radii) 
to the BH event horizon as a function of the Kerr parameter, 
a. Such a solution scheme for accretion onto a non-rotating 
BH has been outlined in Das (2004). We concentrate on ac- 
cretion with sub-Keplerian angular momentum distributions. 
Such weakly rotating flows are exhibited in various physical 
situations, such as detached binary systems fed by accretion 
from OB stellar winds (lUarionov & Sunyaev 1975; Liang & 
Nolan 1984), semi-detached low-mass non-magnetic binaries 
(Bisikalo et al. 1998), and super-massive BHs fed by accre- 
tion from slowly rotating central stellar clusters (Illarionov 
1988; Ho 1999 and references therein). Even for a standard 
Keplerian accretion disk, turbulence may produce such low 
angular momentum flow (e.g., Igumenshchev & Abramowicz 
1999, and references therein). We find that our formalism is 
also valid for accretion with substantial angular momentum 
(as large as Keplerian) for retrograde flows. 

2. FORMALISM 

We consider non-self-gravitating, non-magnetized accretion. 
The gravitational radius is rg = GMbh /c^ and we employ units 
where G = c = Mbh = 1 ■ We use Boyer-Lindquist coordinates 



with signature -+++, and an azimuthally Lorentz boosted or- 
thonormal tetrad basis co-rotating with the accreting fluid. We 
define A to be the specific angular momentum of the flow and 
neglect any gravo-magneto-viscous non-alignment between A 
and a. We focus on the stationary axisymmetric solution of 
the following equations: 

yMcjp.'^O, (py^).^ = 0, (1) 

where 9^^,v'p, and p are the energy momentum tensor, the 
four velocity and the rest mass density of the accreting fluid, 
respectively, and the semicolon denotes the covariant deriva- 
tive. The radial momentum balance condition may be ob- 
tained from (vf^v,j+gfj,,y) Sl^f = 0. Exact solutions of these 
conservation equations require knowledge of the accretion ge- 
ometry and the introduction of a suitable equation of state. We 
concentrate on polytropic accretion for which p = Kp^ , where 
K and 7 are the monotonic and continuous functions of the 
specific entropy density, and the constant adiabatic index of 
the flow, respectively. The specific proper flow enthalpy is 

taken to be /1 = (7- 1) {7- (l + a,) } , where Qs is the poly- 

1 /2 

tropic sound speed defined as flj = {dp/de)g = \l/i(r(r),7) = 
'^zip, Pil)'^ here T{r) is the local flow temperature, e is the 
mass-energy density, and {^'i, ^2} are known functions. The 
subscript S indicates that the derivative is taken at constant 
specific entropy. Hence, a^(r) = jKT(r)/{iJ,mH) = Q^T{r), 
where 9 = [7K/(/im//)]'/^, /i is the mean molecular weight, 
niH is the mass of the hydrogen atom, and k is Boltzmann's 
constant. 

We assume that the disk has a radius dependent local half- 
thickness, H(r), and its central plane coincides with the equa- 
torial plane of the BH. We use the vertically integrated model 
(Matsumoto et al. 1984) in describing black hole accretion 
disks where the equations of motion apply to the equatorial 
plane of the BH, assuming the flow to be in hydrostatic equi- 
librium in the transverse direction. Various dissipative pro- 
cesses are not explicitly taken into account. Even thirty years 
after the discovery of standard accretion disk theory (Shakura 
& Sunyaev 1973), modeling of viscous multi-transonic BH 
accretion, including proper heating and cooling mechanisms, 
is still quite an arduous task, even for Newtonian flows. The 
extremely large radial velocity close to the BH implies r,„ ^ 
r,,,.5, with TiniTyis) the infall (viscous) timescales; hence our as- 
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sumption of inviscid flow is not unjustified, at least out to a 
few or to perhaps tens of rg. One of the most significant ef- 
fects of the introduction of viscosity would be the reduction 
of radial angular momentum. Hence study of the behavior of 
flow variables as functions of A provides some insight about 
the flow behavior for viscous transonic accretion. 

The temporal component of the first part of Eq. (1) leads to 
the conservation of specific flow energy £ along each stream- 
line as, £ = hv, (Anderson 1989). The Kerr metric in the equa- 
torial plane of the black hole may be written as (e.g., Novikov 
& Thome 1973) 



ds^ = gnudx^dx" = -'—^dt^ + ^ {d6-wdtf + ^—dr^+dz^, 
A A 

(2) 

where A = r^-2r + a^, A = + r^a^ + 2m^, and w = 2ar/A. 
The angular velocity f2 is 



(g>0 + ^g») 



4a A(4a^-r^A) 

_r A 

A _4Xa 
T3" f 



(3) 



The normalization relation v^v** = -1, along with the above 
value of fl, gives us the expression 



Vt = 



Ar^A 



(1-m2) {A2-8AaM + AV 



(4) 



where u is the radial three velocity in the co-rotating 
fluid frame. Hence we obtain the conserved specific en- 
ergy (which includes the rest mass energy) to be f = 

(7-i)[7-(i + e2r)]''v,. 

Note that our final calculations will provide the value of u 
and its related quantities rather than those of V(. Hereafter, 
it is implied that any appearance of Vi in any mathematical 
expression is transformed in terms of u using Eq. (4). 

We follow Abramowicz, Lanza & Percival (1997) to derive 
the expression for disk height H{r) in our flow geometry be- 
cause the corresponding equations in their calculations remain 
non-singular on the horizon and can accommodate a thin disk 
geometry as well as quasi- spherical flow. We obtain: 



Hir) = V2r^eT-2 





(5) 

where ip = X^vj - (v, - 1). We define the quasi-invariant en- 
tropy accretion rate, 5, as a quasi-constant multiple of the 

mass accretion rate, 5 = K^Mi„, where dE/dr = along a 
streamUne only for a shock-free non-dissipative flow, and not 
for a flow containing shocks. The mass accretion rate (M) and 
the entropy accretion rate are found to be 



M = 4'KpQT^Mr^ 



7(i-e2rM2) {7- (1+927)} V 



A 



(1-m2)V 



7{7-(l+a2)} 



(6) 



(7) 



We simultaneously solve eqs. (4 - 7) along with the equa- 
tion defining S, to find the dynamical three- velocity gradient 

as 



du 
d? 



(7+1) 


r-1 , 2 v,ax 
A ^ r A'tji 


_ X 

2 




^ 2j/' 


(l-«2) (7+1)(1-«2)k 



(8) 



where 



\ dr dr J 



_ 2 _ 2 A dn 

a-2Xv, — + ^^-^ — ■ (^^^^^^^^^ • 

(9) 

Hereafter, we use the notation [P^] for a set of values of 
{£, A,7} for a fixed Kerr parameter a (including the counter- 
rotating flows where a and A have opposite signs), and [V4] 
for a set of values of {£,X,^,a}. By simultaneously setting 
the numerator and denominator of Eq. (8) to zero we find the 
sonic point conditions: 



as c = 



2ip — u^v,a 



Uc = 



2r(r-l) + 4A 



(10) 



where the subscript c indicates that the quantities are to be 
measured at the critical (sonic) point(s). For a fixed [^3] and 
fl, we substitute the values of Uc and a, |^ in terms of rc (Eq. 10) 
in the expression for £, and obtain a polynomial in rc, the so- 
lution of which provides the location of the sonic point(s) r^. 
To determine the behavior of the solution near the sonic point, 
one needs to evaluate the value of du/dr at that point (the 
'critical velocity gradient' (du/dr)c) by applying L'Hospital's 
rule to Eq. (8). After lengthy algebraic manipulations, we ob- 
tain the following quadratic equation which can be solved to 
obtain (du/dr)c 



du . 



where the coefficients are 
(1 + m2) 2S1S5 



du\ 



(11) 



5i = 



(1-m2)' 7+1' 

a](l-52) 



25\5(, 
74-1 



+ T6, C = -T5\ 



M 1- 



u^VfCF I 2A 

<52 = -r-— , 63 = - + 

lip Vt a 



5a = 62 



2 ^ uvtSo, 



1-m2 



.55 = 



3m2-1 



m(1- 



(7-l-«?)x 



Vt 

I-S2' 



a 

r 

m(7-1- 



1-m2 



Tl = 



73 ■■ 



<--l 2 
aT2X 



2al 

All.! ' 



T2 = 



All) 



^'=A- 



l 2(r-l)2 



2(1-^2)' 

{AX\-a^)ij-v,a^ 

(jijj 
2 VtCT dx 
Alp dr ' 



A2 



T5 = 



7+1 



«?T4 - { (7 - 1 - n + V,a^T3 } 



X 



Idx 
2dr' 



2v,u 



T6 = 



— (7-l-«?)+«?T3 

Vt 



(7+l)(l-M 

Note that all the above quantities are evaluated at r, 



(12) 
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Fig. 1. — A typical flow topology for multi-transonic BH accretion: local 
Mach number plotted against distance from the event horizon of the BH, for 
a = 0.2, 8 = 1.0003, A = 3.1, and 7 = 4/3. The curve ABCD represents 
accretion passing through the outer sonic point B (ro = 1094), while EBI 
represents the self-wind. Flow along GFH passes through the inner sonic 
point F (r, = 3.796) and encompasses a middle sonic point r,,, (r,„ = 14.78) 
denoted by an asterisk. A multi-transonic accretion flow passes through AB, 
then undergoes a shock which causes it to drop to the lower branch of the 
GFH curve. 
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Fig. 2. — Variation of terminal behavior of accretion flow variables, 
M^,ri5,p^, and as functions of BH spin, a. The figure is drawn for a 
shock-free prograde flow passing through the outer sonic point, with E = 
1.00001; the curves (from left to right) correspond to values of A = 2.60, 
2.12, 1.87, 1.55, 1.27, 1.04, 0.829, and 0.609. See text for details. 



3. RESULTS 

For a range of [P4], one obtains three different values of 
for a fixed [V4], among which the largest and the small- 
est values correspond to the X type outer, r^, and inner, r,, 
sonic points, respectively. The O type middle sonic point, r^, 
which is unphysical in the sense that no steady transonic solu- 



tion passes through it, lies between r; and r,,. For > 
and S(r,) < S(ro), one obtains multi-transonic accretion and 
a wind, respectively. If the accretion through r^, can be per- 
turbed in a way so that it produces an amount of entropy ex- 
actly equal to [S(ri)-S(r<,)], supersonic flow through r<, can 
join the subsonic flow through r, by developing a standing 
shock. The exact location of such a shock, as well as the de- 
tails of the post-shock flow, may be obtained by formulating 
and solving the general relativistic Rankine-Hugoniot condi- 
tion (Taub 1948; Thome 1973) for our flow geometry. It is 
worth noting that the idea of shock formation in BH accre- 
tion is contested by some authors (see Narayan et al. 1998 
and references therein for a review). However, the problem of 
not finding shocks in such works perhaps lies in the fact that 
only one (iimer) sonic point close to the BH is explored by 
shock-free ADAF solutions. 

For any [V4], after obtaining {ri,r„,r„], we calculate the 
values of dynamical and acoustic velocities (Eq. 10) and the 
velocity gradient (Eq. 11) for [r,,ro]. These values are eval- 
uated for a specific value of M (= 1.0 MecIcI, in the results 
presented here; however, M could be chosen greater or less 
than the Eddington rate), and we can then compute u{r), as{r), 
the local Mach number, M(r), and T{r). Here we define 
M(r) = u{r)/as{r), but because of our assumption of verti- 
cal equilibrium in an integrated disk model, M(rc) < 1 (Mat- 
sumoto et al. 1984; Chakrabarti 1989). We compute the above 
quantities, the local vertically integrated density, S(r), and 
pressure, W{r) (Matsumoto et al. 1984), and other related dy- 
namical or thermodynamic quantities by integrating the flow 
equations, from and r, down to the immediate vicinity of the 
event horizon, using a fourth order Runge-Kutta method. The 
flow could also be integrated for r > rc to study the subsonic 
part of the flow in detail. Figure 1 shows the integral curves of 
motion for a particular case. Similar topologies are obtained 
for any [Va\ allowing multi-transonic accretion flows. We find 



that r,- 



pro 



< n.r 



for all [7^3] for which valid common solu- 



tions are available for both a and -a; however, for a > 0.2, no 
such common [V^] solutions can be found. 

We define the "terminal" value of any accretion variable 
A at a radial distance rs = (re + 6) as As, where is the event 
horizon, and < i5 <C 1 .0. Here we calculate by integrating 
the flow from I'c down to rs, so is clear that the dependence of 
As on a can be studied at any r. For 5 = 0.01 and Mbh = lOM©, 
Fig. 2 shows the variations of Ms,Ts, ps and ps with BH spin 
for shock free prograde flows passing through rg. 

Note that as a flow through r, does not connect re with in- 
finity (such a flow folds back onto itself, see Fig. 1), solutions 
through r, do not have independent physical existence, and 
can only be accessed if the supersonic flow through rg un- 
dergoes a shock, so that it generates extra entropy, becomes 
subsonic, and produces the physical segment of the solution 
through r, . If a shock forms, the calculation of A5 for shocked 
flows boils down to the calculation of As for flow through 
r,. We find that for flows through r„ the general profiles 
for the variations of A 5 with a remain exactly the same as 
that for flow through r^; only the magnitudes of various As 
change. We obtain Msiri) < Ms(r„) and lTs,ps,ps](i'i) > 
[Ts,ps,ps]{fo), because the shock compresses the flow and 
makes it hotter. Fig. 2 is drawn for 7 = 1.43, but a simi- 
lar figure could be drawn for any value of 7 allowing multi- 
transonic accretion. We find that for both prograde and retro- 
grade flows, {Ts,ps,ps] non-linearly and monotonically cor- 
relate with 7 while anti-correlates with 7. Multi-transonic 
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accretion does not exist for the same [V3] for all values of 
a for prograde flows. Each curve segment in Fig. 2 shows 
the range of a for which a set of {£,X} overlaps to produce 
three sonic points. Specifically, £ is constant in Fig. 2 and 
each curve corresponds to a value of A which allows multi- 
transonicity for that range of a. The allowed overlap range in 
a decreases as a increases. 

For retrograde flows, the situation is different. For many 
combinations of [P3], one obtains a large range of a (some- 
times even covering the complete range of BH spin) for which 
one can have a fixed value of {5, A, 7} producing multi- 
transonicity. Hence the variations of A^'s with a produce con- 
tinuous curves covering almost the whole range -1 < a < 0. 
Therefore multi-transonicity is more common for counter- 
rotating BH accretion. We also found that, unUke prograde 
flows, retrograde flows with values of angular momentum 
even as high as Keplerian can demonstrate multi-transonic be- 
havior. 

From Fig. 2 we see that the the correlations (of tempera- 
ture, density and pressure) and the anti-correlation (of Mach 
number), with BH spin are non-linear but monotonic. Hence, 
rapidly rotating BHs will produce a hotter and denser pro- 
grade flow near the horizon. However, we again found that 
this is not true for counter-rotating flows, where the above- 
mentioned variations are not monotonic; rather they show dis- 
tinct points of inflection at some intermediate value of a. The 
physical reason for such a behavior for retrograde flows is 
not clear to us. If one keeps {S,\, \a\] fixed but decreases A 
(which effectively mimics the introduction of viscosity to the 
system), [Ti,ps,p^\ non-linearly decreases and increases 
for prograde flow, whereas the trends are just the opposite for 
retrograde flow. These trends indicate how our results would 
be modified if our formalism were to be applied to study vis- 
cous transonic flow. Hence the study of A 5 as a function of 
a also opens up the possibility of distinguishing between the 
prograde and retrograde flows if the transonic properties of 
such flows can be measured. 

4. DISCUSSION 
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